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Abstract 

Based on a prior determination of the (ft selfenergy in a nuclear medium we per- 
form a theoretical study of inclusive (ft photoproduction in nuclei, looking at the A 
dependence of the cross sections for different (ft momenta. We find sizeable reductions 
in the nuclear cross sections with respect to the elementary one, using a (ft selfenergy 
which implies a width about six times the free one at normal nuclear density. The cal- 
culations are done to match the set up for an ongoing experiment at SPring8/ Osaka 
which should provide valuable information on the renormalization of the (ft properties 
in nuclei. 



1 Introduction 

The renormalization of the properties of the vector mesons in nuclear matter has re- 
ceived much attention, particularly the p meson, (see for a review). The (ft meson 
has received comparatively less attention, but it turns out to be a well suited tool to 
understand the dynamics of the vector mesons in matter since the changes of some of 
the (ft properties are comparatively larger than those of the p meson, what would make 
its experimental observation in principle easier. The change of the eft properties in the 
medium at finite temperature and/or density has been studied in several approaches like 
effective Lagrangians [111111130 El and QCD sum rules 013 EDI HQ- Particularly, the 
(ft width modification in matter has been subject of study in dropping meson mass sce- 
narios [HllElliniEllinillE! , as a result of collisional broadening through 0-baryon [TJj or 
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0-meson |T8*1 IT§] scattering processes and by modifying the KK decay channel by means 
of in-medium kaon selfenergies [3120] • The majority of these works point at a sizeable 
renormalization of the (fi width and a small mass shift. 

Indeed, in pHITj it was found that the (fi width at normal nuclear matter density, po, was 
about one order of magnitude larger than the free width. Particularly, the results of [JJ 
were based on the use of a K~ selfenergy in nuclear matter which accounted for Pauli 
blocking correction in the intermediate K~N sates in the scattering equation. Further 
studies considering the selfconsistency of the KN in-medium scattering equation were done 
in j2UE2] to determine the K~ selfenergy. Following those results a reanalysis of the (fi 
width in the medium was done in [20], resulting in values around 22 MeV at p = po- A more 
recent evaluation, using the K~ selfenergy of (221 and improving on some approximations 
of [20] was done in [2H1 , resulting in a width of about 28 MeV at p = po an d a small shift 
in the mass of about —6 MeV. 

Several proposals to test the (fi properties in a nuclear medium have been done, for 
instance, regarding the detection of its decay products in A — A and p — A collisions j2UE]- 
Particularly a method to analyse the invariant mass spectra of the emitted charged kaon 
pairs was studied in [2H]. Having into account that the calculations in 001201123 were 
done for a at rest in the nucleus, it has been suggested to determine this large width 
using reactions like ir~p — > (fin in nuclei jZ] and (fi photoproduction in nuclei j2B] , assuming 
the (fi to emerge with a small momentum. The latter proposal is possible by means of the 
elementary 7V — > (fiN reaction with the (fi going in the backward direction and the help 
of the Fermi motion of the nucleons. The small cross section for these particular events 
makes this experiment difficult and in addition it was shown in [2TJ that the distribution 
of the K + K~ pairs from the Coulomb interaction in heavy nuclei removed the changes 
in the (fi width in the K + K~ invariant mass that one might expect from the changes of 
the (fi width in the medium. Another possibility, not previously considered, to investigate 
the changes of the (fi width is to look for (fi loss of flux in (fi nuclear photoproduction. 
The A dependence in this loss is related to the (fi width in the medium. The drawback is 
that in this experiment the largest fraction of the (fi comes out with a momentum of the 
order of 1500MeV/c in the ongoing experiment at SPring8/Osaka [2E1- The combination 
of experiment and theoretical models, which can make the extrapolation of (fi at rest to 
(fi with finite momenta, can thus make the experiment useful to learn about (fi nuclear 
properties. 

With this aim, we present here a theoretical calculation of the A dependence of the 
cross section for (fi photoproduction in nuclei which could test the models of (311711201122] • 
Our work will be based on the models for the (fi selfenergy of Oset-Ramos (OR) [20] and 
Cabrera- Vicente (CV) [23], after a proper extension to account for the finite momentum. 

2 (f> photoproduction cross section 

Let Il(p, p) be the (fi selfenergy in a nuclear medium as a function of its momentum, p, and 
the nuclear density, p. We have 
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^(rj^r/.-^M'O) ( 4 ) 



= V opt = ReV opt + * Im ^opt » (!) 

and hence 

r T n „ imn dp 

- = -Im— ; T = = — , (2) 

2 2u> uj at 

where u is the energy and P is the decay probability, including nuclear quasielastic 
and absorption channels. Hence, we have for the probability of loss of flux per unit length 

dP _ dP _ dP Imn 

~dl ~ Vdt~ ^dl~ p~ ' ^ 
The nuclear cross section for inclusive photoproduction will then be 

da A [ , 3 -, , sda 

ln = J drp{r) dn 

where 4§ and are the elementary and nuclear differential cross section and f' = f + 
with r the production point inside the nucleus. The exponential factor in Eq. (j3J) 
represents the survival probability of the meson in its way out of the nucleus. If the 
did not get absorbed inside the nucleus then we would get the typical result for an 
electromagnetic reaction -Sf = Aj^. Eq. (0J) relies in the eikonal approximation which is 
accurate for the large momentum involved in the process. In practice one might expect 
small corrections to the formula, even if the did not decay, from two sources: 

i) Distortion of the trajectory because of the real part of the potential 

ii) Change of direction and energy of the m quasielastic collisions <pN — > <p'N'. 

The first effect should be negligible since the real part of the potential is so weak that 
only modifies the mass of the in about 6MeV jSlEBj- The second effect should be very 
weak too, account taken of the lack of direct coupling <j)NN because of the OZI rule. On 
the other hand, the effect from this source would simply lead to a slight change in the 
direction but not the disappearance of the 0. This simply means that collecting the in a 
narrow cone along the forward direction, where practically all the go, both because of the 
kinematics of the lab variables and the extremely forward direction of the 7iV — > <pN cross 
section [21], would guarantee that the undergoing quasielastic collisions in the nucleus 
are accounted for in the nuclear cross section. In order to adjust to this experiment we 
should not remove theoretically the events in which there are quasielastic collisions. This 
will be done by including in Imll only the absorption events. 

The integral in Eq. (J3J does not depend on the direction of the momentum and thus 
we have 

P out = ^L = LJd>r P (r)el~*i™^ (5) 

and this is the magnitude which we would evaluate as a function of p and A, which can be 
interpreted as the probability for a to go out of the nucleus. The density profiles, p(r), 
for the different nuclei have been taken from [3*U] . 
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3 Evaluation of the <\> selfenergy at finite momentum 



For the evaluation of Imll in nuclear matter we are going to use the OR and CV models 
which are based on chiral SU (3) dynamics considering the K and K in-medium properties. 
In these works the evaluations were done for a at rest. Since in the present work we 
are dealing with mesons with momentum up to 2 GeV, an extension of the models to 
finite momenta is mandatory. In the following subsections we summarize these two models 
stressing the modifications to consider the finite momentum. 

3.1 Extension of the Oset-Ramos model to finite momentum 

The OR accounts for —>■ K + K~ and — > K°K° decay diagrams where the kaons 
are renormalized in the medium due to s-wave and p-wave interactions, Figs. Q and El 
respectively. 




a) b) 



Figure 1: Diagrams contributing to the selfenergy coming from the s-wave KN interac- 
tion. Y represents A, S or 5. 




^ p-q 

K 

Figure 2: Diagrams contributing to the selfenergy coming from the p-wave KN interac- 
tion. Y represents A, S or S*(1385). 

To the diagrams of Fig.Elone must add the vertex corrections shown in Fig.|3]for consistency 
with gauge invariance. 

In the evaluation of the previous diagrams, the K + , K°, selfenergy is taken as in |3*T| 
following the tp approximation 
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Figure 3: Diagrams accounting for vertex corrections. 



U K+ K0 = 0.13 mt-f- MeV . (6) 
Po 

Concerning the K, a large contribution to the <fi selfenergy comes from the p-wave 
and vertex correction diagrams. On the other hand, since the K~ p-wave selfenergy is 
proportional to q 2 , it is likely to be more affected inside the loops by the finiteness of 
the <fi momentum. Hence we evaluate explicitly the selfenergy from the p-wave K~ 
selfenergy diagrams at finite cf) momentum and take the same part coming from the K~ 
s-wave selfenergy as obtained for zero momentum. This latter part produces a contribution 
to the imaginary part of the <fi selfenergy that can be parametrized as 

Imn (s) = -7.Q2p/p uj(p) MeV (7) 

We now evaluate the contribution to Imll(p, p) of the K p-wave selfenergy diagrams, 
depicted in Fig. [21 

In Fig. K represents K + or K°, K represents K~ or K° and Y the hyperons A, £ 
and S*(1385). The <j) selfenergy coming from these diagrams at finite momentum can be 
obtained as 

-fllW(p,p) = 2 J ^JD(q)iD(q)(-i)uf(q,p)iD(p-q lP )^V^ R V^ KR = 
= H J ^D(qfrt?\q,p)D(p-q,p^ 

where D(q) = (q 2 — m 2 K + ie)^ 1 is the K~, K°, propagator, D{p — q) — ((p — q) 2 — m 2 K — 
n^+) -1 is the K + , K°, propagator with the selfenergy of Eq. ©. In the evaluation of 
Eq. (jHJ) we have used that V^kk — ~W<j) e ^{PK ~ Pk)-> 9<t> = 4.57. The factor 2 at 
the beginning of the equation is due to the possibility of having K~K + or K°K° in the 
intermediate states. 

To obtain Imll we can apply the following Cutkosky rules for the cut represented by 
the dotted line of Fig. 

U(k) -> 2iImn(A;) 
D(k) -> 2iQ(k°)lmD(k) 



(q-pf 
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which leads, after performing the q° integration, to 



Imn (p) (p, p) 



(2tt) j 



m 



K 



to. 



p-q 



(q ■ pf 



- q 



9(g°) 



(9) 



with to. 



p-q 



P ~ 



Jp - q and 9 is the step function. The 

K p-wave selfenergy, ImlT^g), is evaluated in 22J and is a function of the KNY vertices 
and the Lindhard functions for hyperons. In Eq. Q we also add, as in Ref. j?j, a form 
factor for each kaon-baryon vertex of dipole type, [A 2 /(A 2 — q 2 )} 2 , with A = 1.05 GeV. 

Finally, the evaluation of the vertex corrections of Fig. |3] is done in an analogous way. 
In this case we need the <pKNY vertex functions given by 



V<f>KNY 



9<t> V KNY & ■ e 



(0). 



(0) 



Y 



A, E 



(10) 



where Vkny an d Vrnt,* are coefficients given in [20J. The resulting expression for the 
imaginary part of the (f) selfenergy coming from the vertex corrections is the same as 
Eq. © but doing the following substitution 



2^2 



D(q) 2 q 



(q-pf 



mZ 



D(q): 



^2 _ (p-q)(p- q) 

q i 

ml 



+ 1 + 




3.2 Extension of the Cabrera- Vicente model to finite momentum 



The (f> selfenergy from KK loop diagrams, calculated in [71122] for a at rest in a nuclear 
medium, is given by 



d A q 



(p ■ qf 2 

- q 



mZ 



D K (P° -q°,P~ T, p)D R (q°, q; p) (12) 



for a meson with a momentum p. 

A spectral representation [7j which sums up to all orders the insertion of irreducible 
kaon selfenergy terms is used in the kaon propagator. The imaginary part of IKf , which 
is of interest in the present work, can be written as 



1 o4 



Imn£*(p°,p;p) = -— sj- / dqq 
(p-q) 2 



du— — 

1 ^p-q 



q°=P°-UJp-q 



Sr{p° -Up- q ,q;p)®(p° -u p - q ) , (13) 
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Figure 4: Imaginary part of the in medium selfenergy, without the inclusion of the free 
part, as a function of the momentum of the and the density, for the finite momentum 
modification of model [2U] (dashed line) and 23; (solid line). 



The main differences with respect to the approach described in section 13.11 are the 
following: 

• Eq. (fT5"j) considers higher order effects in the insertion of the kaon selfenergies as 
compared to Eqs. (J7J) and Q. Moreover, the calculation at finite momentum is considered 
for both the s- and p-wave K selfenergy contributions. Note also that Eq. (|T3j) includes 
the free width into the KK channel, which will be subtracted to keep only the nuclear 
effects in the selfenergy. 

• We use the p-wave K selfenergy given in j2H] which, starting from the result in [2Ti] . 
takes into account an improvement of the relativistic recoil corrections of the KNY vertices. 
In addition, these vertices carry a static form factor of the form [A 2 /(A 2 + q 2 )} 2 , with 



• The vertex correction diagrams considered in section |3~T1 are also calculated here at 
finite momentum, with the prescription of using the fully dressed kaon propagators in 
the medium. 




(14) 



A 



1.05 GeV. 



4 Results and discussion 



In Fig. 0] we show the imaginary part of the selfenergy as a function of the momentum 
for different nuclear densities and for the two models described in section 13.11 (dashed line) 




Figure 5: a^/^Aa) as a function of the nuclear mass number for two different momentum 
of the calculated with the two models described in the text. 



and 13.21 f solid line). The contribution to Imfl coming from the free decay into KK, non 
density dependent, has been subtracted from the full Imll since the KK coming from the 
free decay would be detected and counted as a event, hence it does not contribute to 
the loss of flux required in the argument of the exponential in Eq. (|5|) . Despite the visible 
differences in Fig. H] between the two models at zero momentum, the trend of the plot 
as the momentum increases is very similar. This initial discrepancy at zero momentum is 
expectable due to the differences in the treatment of the s-wave contribution, the relativistic 
recoil corrections and the fact that the model of section 13.21 goes beyond the first order in 
density, implicitly considered in the model of section 13.11 The differences between the two 
models in Fig. IH however, are indicative of the intrinsic theoretical uncertainties. 

In Fig. El we show the ratio <ja/{A<j), which represents the probability of one photo- 
produced to go out the nucleus, calculated theoretically as a function of A and for two 
different momenta. Again, a comparison between the two models is shown. We observe 
that for p = 1400 MeV and heavy nuclei this ratio can be of the order of 0.65, indicating 
a clear loss of flux which should be identifiable experimentally. 

In Fig. Elwe show instead the results for <ja/(A<j) as a function of the momentum for 
two nuclei and comparing the two theoretical models. In this figure we observe that the 
amount of flux lost is larger for smaller momentum. This is logical to expect since the 
probability of decay per unit length is — Imll/p, which is larger for small momenta because 
of the factor 1/p. This occurs in spite of the fact that —Imll decreases with p, as can be 
seen in Fig. |U because this p dependence is weaker than that of 1/p. The set up of the 
experiment at SPring8/Osaka |2Hj is such that it produces momenta around 1500 MeV, 
using photons from 1.4 to 2.4 GeV, coming from production in the forward direction 
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Figure 6: o^ji^Aa) as a function of the momentum of the 
lated with the two models described in the text. 



for two different nuclei calcu- 



in the CM frame. A possibility to extract further information is to use photons of lower 
energy, around, 1.6 GeV, not far from threshold which leads to <fi of around 100 MeV in the 
CM frame and around 1000 MeV in the lab frame. As one can see in Fig. El in a nucleus 
like 64 Cw the depletion factor goes from 0.8 at p = 2000 MeV to 0.7 at p = 1000 MeV. 
This latter factor would be 0.57 for A ~ 240. 

So far the only nuclear effects considered have been the ones related to the absorption of 
the meson. At this point it is worth mentioning that other nuclear effects regarding the 
production mechanism may lead to a further loss of flux or change in the distribution. 
These other nuclear effects are mainly the Pauli blocking of the final nucleon and the Fermi 
motion of the initial one. The first one may lead to a reduction of the flux in comparison 
to the free case because a certain amount of events are forbidden due to the Pauli blocking 
of the final nucleon. On the other hand the initial Fermi motion can distort the distribution 
of the final mesons. In order to estimate the possible flux reduction due to these sources, 
we have included in the integrand of Eq. (jSJ) a factor G(Q,p) which considers a Fermi 
average of these effects [32"] : 

G(Q,p) = 1 - 9(2 - Q) (l - ~Q + ^Q 3 ) (15) 

where Q = \Q\/kp with Q the momentum transfer of the nucleon and kp = (|7r 2 p(r)) 1 / 3 
is the Fermi momentum of the nucleons. 

In Fig. we show aj\,/(Aa), with the CV model, as a function of the momentum 
for m Cu. In solid line we show the result without considering the effect of Eq. (|15|) . (i.e., 
the same as in Fig. |BJ). The other lines represent the results considering the Pauli effect 
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Figure 7: Effect of the Pauli-blocking correction in the ca/{Ao) as a function of the 
momentum of the 6 for 64 Cw calculated with the CV model. 



estimation, using Eq. (fTKJ) . Note that the Q dependence of this G factor introduces another 
kinematical variable in the P ou t observable. We have chosen this extra variable to be the 
zenithal angle of the meson in the Lab frame {0^ ab ). The dashed and dot-dashed curves 
represent the results considering the Pauli effect for 8^ ab = 0° and 5° respectively. We 
can see that significant effects are obtained for large momenta, what is expectable since 
large momenta imply small final nucleon momenta, which are strongly Pauli blocked. 
The Pauli effect also decreases with the meson angle because the momentum transfer 
strongly increases with this angle, therefore the largest effect is obtained for the forward 
direction. Since the photoproduction is strongly forward peaked, one can expect P out to 
be actually somewhere between the dashed and dot-dashed lines in a real experiment. 

In Fig. Owe show the A dependence of P ou t, with the CV model, without including the 
Pauli correction (solid line) and including it for 8^ ab = 0° and 5° (dashed and dot-dashed 
lines respectively), all for = 2000 MeV where the effect is maximum. We observe that 
the loss of flux due to the Pauli correction is nearly constant in a wide range of the mass 
number, resulting in a reduction of around 0.1 in P out for the forward direction. Nonetheless 
this effect is relatively smaller for heavier nuclei compared to the absorption. 

The cross sections calculated in this work are inclusive, summing over all possible final 
nuclear excited states. The coherent cross section, where the final nucleus is the original 
one, is not included in the calculation. Its evaluation requires a complete knowledge of 
the spin and isospin dependence of the elementary 7V — > 0iV amplitude and a different 
treatment than the one done here. The coherent cross section involves the square of the 
product of the mass number and the nuclear form factor times the square of the spin- 
isospin averaged amplitude. The coherent and incoherent cross sections can in principle be 
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Figure 8: Effect of the Pauli-blocking correction in the a a/ (Aa) as a function of the nuclear 
mass number for = 2000 MeV calculated with the CV model. 

separated experimentally, as done in the case of pions around the delta energy region [33] . 
However, this might be difficult at higher energies where one usually has a poorer energy 
resolution. An alternative to this separation, concerning the present work, is to look for </> 
production at finite angles where the momentum transfer to the nucleus is large and the 
coherent process can be neglected. The incoherent cross section sums the square of the spin- 
flip and non spin-flip amplitudes, while in the coherent part only the spin independent part 
contributes. Thus, the ratio \A F{Q)\ 2 /(AG(Q)), where F(Q) is the nuclear form factor 
and G(Q) the Pauli blocking factor of Eq. f(15j) . is an upper bound for the ratio of coherent 
to incoherent nuclear cross sections. We have studied this ratio as a function of the photon 
energy, the </> meson angle and the mass number. Our results can be summarized as follows: 
1) We find that below _E 7 = 2000 MeV the ratio of coherent to incoherent cross section is 
smaller than ten percent for all angles and nuclei beyond A ~ 12. 2) The ratio increases 
rapidly with the photon energy, since the momentum transfer decreases. For instance at 
E 1 = 2400 MeV, for l2 C and forward angles, the ratio can be of the order of 50 percent, 
but for heavier nuclei the ratio is smaller. For instance, for 21 Al this ratio is smaller than 
10 percent. We also obtain that beyond 4 degrees the ratio is only a few percent and can be 
safely neglected. In practical terms, from the experimental point of view, we can say that 
in nuclei around l6 or heavier and photon energies smaller than 2000 MeV, the coherent 
contribution is negligible. This discussion may serve to select an experimental setup in 
which the coherent contribution can be neglected, hence facilitating the interpretation of 
the data. 
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5 Conclusions 



We have shown that using present models for the (ft selfenergy in a nuclear medium, con- 
veniently extrapolated to finite (ft momenta, it is possible to evaluate the survival rate of (ft 
produced in nuclear (ft photoproduction, and how this survival rate is tied to the (ft width 
in the medium and the momentum of the (ft. The survival rates for (ft coming from photons 
in the range of 1.6 to 2.4 GeV are of the order of 0.7, a significant deviation from unity, 
which are measurable experimentally. We have shown the A dependence of the expected 
<ja/{A<j) ratio as well as its dependence on the (ft momentum. Comparison of the results 
with experimental numbers of the incoming experiments would help determine the accu- 
racy of the models used. These models could then be used to extrapolate results at other (ft 
momenta and one could then get a fair idea of how the (ft properties are modified in nuclear 
matter for (ft at any finite momenta below those studied in the present work. 
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